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COMPUTATION OF THE TWO REGULAR 
SUPER-EXPONENTIALS TO BASE EXP(1/E) 

HENRYK TRAPPMANN AND DMITRII KOUZNETSOV 



Abstract. The two regular super-exponentials to base exp(l/e) are con- 
structed. An efficient algorithm for the evaluation of these super-exponentials 
and their inverse functions is suggested and compared to the already published 
results. 



1. Introduction 

We call a holomorphic function F a superfunction [1(1 of some base function / 
if it is a solution of the equation 

(1) F{z + l) = f{F{z)) 

In the case / — exp^, i.e. for the exponential base function f{z) = b', we call F 
super- exponential to base b. In addition, if the super-exponential F satisfies the 
equation 

(2) F{0) = 1 

we call it tetrational; for integer values of the argument z, equation ([T]) and ([2| 
implies F to be the z times application of the exponential expf, to unity 

(3) =expJ expJ...expf,(l).. 



Conversely a function A is called Abel function of some base function / if it 
satisfies 

(4) A{f{z))=A{z) + l. 

For f{z) = b^ we call A super-logarithm to base b. The inverse of a super-exponential 
is a super-logarithm. (In some ranges of values of z, the relations F{A{z)) = z and 
A{F{z)) = z hold.) 

We have constructed super-exponentials and efficient algorithms of their numer- 
ical evaluation in [7] for 1 < & < e^/°, and in [9J for 6 > e^/°. However, neither 
method used in these publications is applicable to base 6 = e^/°. Especially this 
case is analyzed by Walker in [18]; he evaluates the two Abel functions at several 
points in the complex plane. Here we show that his constructions are equal to 
the two regular Abel functions {regular in the sense of Szekeres [TS]) and suggest 
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Figure 1. Super-exponentials to base b = exp(l/e) versus real 
argument; the asymptotic y = e; circles represent the data from 



a faster/more precise alternative algorithm (which goes back to Ecalle) that allows 
to plot the complex maps in real time. 

The two super-exponentials Fi and Fj, along the real axis are shown in figure [l] 
The circles represent the data from tables 1 and 3 by [l8j. The behavior of these 
functions and their inverses in the complex plane is shown in figure [2j 

As in 9J, the subscript of the super-exponential (here 1 or 3) indicates the value 
at of the chosen representative of the class of all super-exponentials obtained by 
argument shift F(z + c), c e C . We often identify this whole class as one super- 
exponential. We consider two classes of super-exponentials represented by Fi with 
Fi(0) — 1 and by F3 with ^3(0) = 3, respectively. According to the definition, Fi 
is a tetrational. For the other (above unbounded) super-exponential, the smallest 
integer from the range of values along the real axis is chosen as value at zero. 

2. Four methods of calculating the regular iteration with 

multiplier 1 

In the theory of regular iteration (see e.g. [15 or T2]) there are several algorithms 
available to compute the regular fractional/continuous iteration and the Abel func- 
tion of an analytic function at the fixed point 0. Functions h with multiplier 1, e.g. 
h!{0) — 1, are treated differently from functions h with |/i'(0)| 7^ 0, 1. 

In our case we have the base function f{z) = e^/° with fixed point e and /'(e) = 1. 
As the whole theory of regular iteration assumes the fixed point to be at 0, we 
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Table 1. Computing the Abel function with Levy's formula. 
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n 


Vn 


n 




Vn 


n 




Vn 


n 




Vn 


100 


-1.4560 


1,000 


-1 


425788 


10,000 


-1 


4226982 


100,000 


-1 


42241848 


101 


-1.4557 


1,001 


-1 


425785 


10,001 


-1 


4226982 


100,001 


-1 


42241893 


102 


-1.4553 


1,002 


-1 


425781 


10,002 


-1 


4226981 


100,002 


-1 


42241823 


103 


-1.4550 


1,003 


-1 


425778 


10,003 


-1 


4226981 


100,003 


-1 


42241951 


104 


-1.4547 


1,004 


-1 


425775 


10,004 


-1 


4226981 


100,004 


-1 


42241880 


105 


-1.4544 


1,005 


-1 


425771 


10,005 


-1 


4226980 


100,005 


-1 


42241891 


106 


-1.4541 


1,006 


-1 


425768 


10,006 


-1 


4226980 


100,006 


-1 


42241937 


107 


-1.4538 


1,007 


-1 


425764 


10,007 


-1 


4226980 


100,007 


-1 


42241983 


108 


-1.4535 


1,008 


-1 


425761 


10,008 


-1 


4226979 


100,008 


-1 


42241913 


109 


-1.4533 


1,009 


-1 


425758 


10,009 


-1 


4226979 


100,009 


-1 


42241958 



move the fixed point to via a conjugation with the linear transformation r: Let 
t(z) = e (z + 1) then t~^{z) = z/e — 1 and 

T-^ofoT{z) = c^ - 1 =: h{z) 

f = T O h O 

The regular iterates at the fixed point e are then given by = r o /il*! o t~^, 
where the regular iterates of h can be obtained in one of the later described ways. 
The regular Abel function a of h (up to an additive constant determined by u) is 
defined by the inverse of cru(t) = We call this (T„ the regular superfunction 

of h with fTi, (0) = u. 

(5) fT„(i) a„:=(7-i /jl*! (z) = + a„(z)) 

The regular Abel function A„ of / at e with Au{u) = and the regular super- 
function Fu of / at e with Fu{0) = u can be obtained by 

(6) Au = ar-i{u}°T^^ Fu^Toa^^i(^y 

The classic limit formula of Levy [T3] (see also Kuczma pl2j theorem 3.5.6) for 
the regular Abel functions of h with multiplier 1 is: 

(7) au{z) = hm ' ' ^ ' 



n^oo /l[«+l](u) - /?,[«] (u) 

One can verify that in our case of h{z) =c^ — 1, Levy's formula converges just too 
slowly; it is difficult to reach sufficient precision to make any camera-ready plot of 
the Abel function. In table [T] we display 

_ /[«](_!)_ /M(l) 

/[n+l](l)_/N(l) ^^1^ 

There is another interesting possibility to compute the regular superfunction, 
which we call here Newton limit formula (probably first mentioned by Ecalle in [5]) 
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Table 2. Computing the Abel function with Fatou's formula. 



n 


an 


n 


On 


n 


a„ 


1,000 


-1.4224939 


10,000 


-1.422367740 


100,000 


-1.42235507550 


1,001 


-1.4224938 


10,001 


-1.422367738 


100,001 


-1.42235507549 


1,002 


-1.4224936 


10,002 


-1.422367737 


100,002 


-1.42235507548 


1,003 


-1.4224935 


10,003 


-1.422367736 


100,003 


-1.42235507546 


1,004 


-1.4224934 


10,004 


-1.422367734 


100,004 


-1.42235507545 


1,005 


-1.4224932 


10,005 


-1.422367733 


100,005 


-1.42235507543 



because of its similarity to the Newton binomial series of = (1 + (a; — 1))*: 



(10) F„(t) = /[*](.) = E E (^) (-i)"-"7""i("). 

n=0 ^ ^ m=0 ^ ^ 

However also this method has a depressing slow convergence, moreover we need a 
bigger internal precision caused by the involved summation. For example for 1000 
summands and 2000 bits precision with u = 1 and t = —1.4223536677333 we get 
Fu{t) ~ —0.9875 while we would expect a value very close to —1 (see table[T]). 

Another formula to compute an Abel function of e^ — 1 is given in Walker's text 
jl8j . He computes an Abel function gi of a{z) = 1 — e~^ and an Abel function 52 
of h with a formula which goes back to Fatou [4]: 

1 2 

(11) gi(z)= lim --log(n)+ - n, z < 0. 

n-i-oo ^5 ai"i (z) 

We derive the corresponding Abel function of h by knowing that a(z) — —h{~z). 

1 2 

(12) a\^\z) ^ gi{-z) = Ihn^-- login) - -^-^-n, z<0 

1 2 

(13) «w (z) = 92iz) ^ lirn^ -- log(n) - ^[_„]^.^-^ + n, z>0 

The convergence of this formula is better than that of Levy but still rather slow 
(which Walker notices too and that's why he introduces a slightly accelerated ver- 
sion which we omit here for brevity). To have an impression of the convergence of 
Fatou's/Walker's formula, we display 

2 2 

y^'- /jH(_e-i-l) + /iN(_l) 

^ aW(r-i(-l)) - a{^\r-\0)) 1 = ^i(-l) 

in table [21 

Before we give the fourth method and showing that Walker's formula is equivalent 
to it, we start with some formal background about regular iteration. 

Definition 1 (regular iteration). For every formal powerseries 

00 

(14) h{z) = z + ^ hnz'\ m>2,h^^0 
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and each t e C there is exactly one formal powerseries h'^*^{z) ~ z + J2'^=m 

such that h^*\n = t ■ hm and h'^*^ o h — h o h'^*^ . We call ft.'*' the regular iteration of 

h. It satisfies h^^^ — h and — ft.I'^l o /il*! and is given by the formula: 



N-l , 

t \ \ ^ / n 



n—Q ^ m— 

Af-1 



— 1 — m 



N 



(16) -E(-i 



m=0 



N-l-i 



N 



where (161 can already be found as formula (2.19) in 

The formal powerseries /i'*! is not necessarily convergent even if h is. We call 
a function which has the powerseries as asymptotic expansion at a regular 
iteration of h. 

If z I— >■ (z) is an analytic function in some domain, we call the function 
<7{t) — (zq) a regular superfunction of h for any zq in the domain, and we call its 
inverse a regular Abel function of h. Usually we identify Abel functions that only 
differ by a constant and we identify superfunctions that are translations of each 
other (i.e. x i— > cr{x + c) is identified with a). 



Definition 2. Let h he a. formal powerseries of the form h(x) — x + J2T=m ^kX^ i 
hm 7^ 0. Its iterative logarithm is the unique formal powerseries j of form j{x) ~ 
YlV=m^kX^ with jm — h„i that satisfies the Julia equation 

(17) joh^h'- j. 

One obtains the Julia equation when differentiating the Abel equation and then 
substituting a' — l/j. For reference we give the first few coefficients of the iterative 
logarithm j of h{x) = e^ — 1: 

V N 1 2 1 3 1 4 1 5 11 6 1 7 

i(x) = -X x'^ H X H x^ x' + . . . 

■'^ ' 2 12 48 180 8640 6720 

From this iterative logarithm one can get a description of the regular Abel function 

= = + - + 270" + 2^^' - 10^ 

If we integrate this to get a the term x^^ becomes log \x\ for real x and log(±a;) for 
complex values of x; the choice of the sign determines the branch of the resulting 
function (that unavoidably has the cutline). This gives the expansion 



s{x) v{x) 

This formula can not be used to get arbitrary precision, because v{x) is not con- 
vergent as we show now; however suitable truncation of the divergent series can be 
used to obtain a certain precision. 

Preliminary 1 (Baker 1958 [T] Satz 17). The regular iteration /i'*! of h{x) = e^ — 1 
has non-zero convergence radius exactly if t is an integer. 
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Preliminary 2 (Ecalle 1975 3J ). Let h he a formal powerseries with multiplier 
1. Its regular iteration powerseries /i^*! has a positive radius of convergence for all 
t (z C if and only if its iterative logarithm has a positive radius of convergence. 



Theorem 3. The formal powerseries v[z) in (18) has convergence radius. 

Proof. Suppose that v{z) has non-zero radius of convergence. Then also v'{z) = 
a'{z) — s'{z) has non-zero radius of convergence. Then z'^v{z) has non-zero radius 
of convergence and then z^s'(z) + z'^v(z) is a powerseries with non-zero radius of 
convergence with non-zero zcroth coefficient. Then 

a'{z) z^s'{z) + z^v{z) 

has non-zero radius of convergence. Then by theorem [2] the regular iteration /i^*! of 
h{x) = e^ — 1 has non-zero radius of convergence for all t. This is in contradiction 
to theorem [1] Hence, the series diverges. □ 



Nonetheless we can use the formula ( 18 1 in a different way to calculate the regular 
Abel function of h(x) — e^— 1. If we truncate v to N summands, denoted by vn and 
aN '■= s+vn, then Ecalle showed (see [2], P- 78 ff., 95 ff.) that there are 2(m— 1) 
different regular Abel functions a'^^\ 1 < J < 2(m — 1) (a^^^ is defined on the j-th 
petal — each petal touching the fixed point — of the now called Leau-Fatou 
flower see [2]) given by: 

(19) a(j')(z) = lim a^^(/l(-i)'^'"l(z)) - (-1)^+171, NeN 

n— f 00 

where a^''(z) is aN{z) with the logarithmic term logjzl (in (I81) replaced by 

log l^ze"'-^"'^-' "-1 •''^ and 9q is the unique number in the interval {—i^^^, T^T^] such 

that km = |/i™|e-'('"-i)««. 

This apphed to h{z) = e^ — 1 where m = 2, km = 1/2, 9o = 0, we get the two 
regular Abel functions 

(20) a(i)(z) = jim i log(-;iN(^)) _ _2_ + v^{h^^Hz)) - n, 5i(z) < 



n— J-oo 



(21) a<^^Hz)^ lha\log{h^--^iz))--j^ +v^{h^--^{z))+n, 5R(z) > 



< e 



(22) yl(i)(z) = a(i)(^-l),3fJ(z) 

(23) (z) = - 1) , ^{z) > e 



From these 4 methods we know that Levy's formula Q, the Newton formula ^ 
and Ecalle's method (20) and (21) calculate the ref^wZar iteration/ Abel function. 
We show in the last part of this section that Walker's formula is equal to Ecalle's 
formula and hence also computes the regular Abel function. 

The application of theorem 1.3.5 in [T^ gives the following: 

Preliminary 4 (Thron 1960 [TB] Theorem 3.1.). Let h be analytic at with pow- 
erseries expansion of the following form 

h{x) ^x + /i,„a;™ -|- + . . . , /i™ < 0, m > 2 
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then 

n— ^CJO 

Now, about the functions and a'^ constructed by Walker, we have the 
following theorem: 



Theorem 5. Functions and a'^ given in (12) and (13) are the two regular 
Abel functions of x i— > e^ — 1. 

Proof. We show that the difference of Walker's and Ecalle's limit formulas is a 
constant. The differences are: 

Si{z) = lim i log(-/i["l(z)) + l log(n) = lim ^ log (~Tih^"\z)) , z < 

n— yoo 3 3 ri— >oo 3 V / 

62{z) = lim - log(/i[-"l(z)) + - log(n) = lim ^ log f n/il""! (z)) , z > 

ri— >oo 3 3 ri— >oo 3 V / 

As X i-> ~h(—x) and x i— >■ for h{x) = e^ — 1 is of the form required by 

preliminary UJ with to — 1 = 1 we see that each of n/i["l(z) and nh^~"\z) converges 
to a constant independent on z. □ 

3. A NEW EXPANSION OF THE SUPER-EXPONENTIALS 

This section describes an evaluation of the two super-exponentials to base b = 
exp(l/e). The given expansion is fast and precise; it allows to plot the complex 
maps of these functions in real time. These maps are shown in figure [2] 

The base function h{z) = expj(z) = exp(z/e) has the only fixed point z = e. 
The super-exponential is expected to approach this point asymptotically. Consider 
the expansion of the super-exponential / in the following form: 



(24) F{z 



2 ^ P,„(-ln(±z)) /|ln(z)| 
zi (3z)™ ^ z^+i 



m+l 



where 

m 

(25) Pm{t) = Y,Cn,mr 

The substitution of ([T]) into equation 

(26) F(z-f 1) = exp(F(z)/e) 

and the asymptotic analysis with small parameter |l/z| determines the coefficients 
c in the polynomials ( 25 ) . In particular, 

t 

t^ + t+l/2 

^2^ 2^ 10 

4 13 o 45 . 53 67 

r H — r H — r h — 1-\ — 
3 6 10 60 

. 77 4 101 o 83 2 653 2701 

f H 1^ H 1^ H r H 1 H 

12 6 4 60 1680 



(27) 


Piit) 


(28) 


P2{t) 


(29) 




(30) 


Piit) 


(31) 


P^it) 
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-8 -6 -4 -2 2 4 6 



10 12 14 16 18 20 22 24 26 28 x 



Figure 2. Map of f — Fi{z), top, and f — F3{z), bottom, in the 
plane z = x + iy. Levels p = 3fi(/) = const and q = — const are 
shown; thick lines correspond to the integer values. 



BASE EXP(1/E) 



9 



The evaluation with 9 polynomials P gives an approximation of /(z) with 15 deci- 
mal digits at 3fi(z) > 4. For small values of z, the iterations of formula 

(32) F(z) = ln(i^(z + l)) e 

can be used. With complex<double> precision, the resulting approximation re- 
turns of order of 14 correct decimal digits in the whole complex plane, except the 
singularities. 

For the tetrational we choose the negative sign inside the logarithm t = — ln(— z). 
Then 

(33) Fi{z)^F{z + xi) 

where xi ~ 2.798248154231454 is the solution of the equation f{xi) = 1. For real 
values of the argument, this function is shown at the bottom of figure [T] The 
complex map of this function is shown at the top of figure [2] 



The same expressions (31 1 can be used also for the above unbounded super- 
exponential with t — — ln(z). Then, the expression 

(34) F(z) = exp(i^(z-l)/e) 

allows the evaluation of the above unbounded super-exponential at small |z|. The 
specific super-exponential can be expressed as 

(35) F:i{z)^F{z + x^) 

where 2:3 sa -20.28740458994004 is solution of equation ^(zs) = 3. Function F3 
by (35),([34|,([T|) is shown at the top of figure [I] for real argument and in the bottom 



picture of figure [2] for the complex values of its argument. 

Due to the leading term in the asymptotic representation, at large values of |z| 
(except the vicinity of the real axis), both functions ^3(z) and i^i(z) behave similar 
to the function z 1/z. 

F3 is entire, and shows fast growth along the real axis. At large values of |z|, the 
function -F'i(z) approaches value e. Function ^3(z) approaches value e for |z| 00 
except on the positive direction of the real axis; in this direction, this function 
shows "faster than any exponential" growth. 

4. Numerics and behavior of the two super-logarithms 



For the evaluation of the super-logarithms, the expression ( 18 1 is used. The 
truncation of the series v keeping the term of the 15th power was used to build-up 
an approximation of A^^^ that returns at least 15 decimal digits for |z/e — 1| < i. 

(36) A(^\z) = a(-C)«^^ + ^ + Ec„C 

where C = (e — z)/e. For larger values, the representation 

(37) A(i)(z) = A(i)(exp(z/e)) +1 

is iteratively used. This allows to extend the approximation to a wide domain 
keeping of order of 14 correct decimal digits. Then 

(38) Ai{z) = A^^'> (z) ~ A^^\l) « ^(i)(z)- 3.029297214418 

is the regular Abel function with the additive constant chosen such that ^i(l)=0. 
This function is shown in the top of figure [3j 
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-8 -6 -4 -2 2 4 6 8 10 12 14 16 18 20 22 24 26 28 X 

Figure 3. The Abel functions f — Ai{z) , top, and J — Az{z), 
bottom, in the same notations as in figure [2j 



BASE EXP(1/E) 



11 



The function Ai is periodic; the period is Ti = 27rei « 17.079468445347134131 i. 
For real values z>e, the representation diverges, indicating a natural way to place 
the cut of the range of holomorphism. In vicinity of this cut, the Abel- function 
shows complicated, fractal-like behavior: the self-similar structures reproduce long 
the range with high density of levels of constant real or imaginary part of Ai. 

At the left hand side of the picture, Ai{z) approaches its asymptotic value —2 as 
3?(z) —oo . Along the strips in vicinity ^{z) ~ 5(Ti -I- 7ri(27i-|-l)), nG integers, 
as 3?(z) — )• +00, this function approaches its another limiting value: Ai{z) — > —3. 
The transfer from the asymptotic value —2 to the asymptotic value —3 corresponds 
to the transition from the singularity at z — —2 to the singularity z — —3 of function 
Fi{z) in the top picture of figure [2] 

The second super-logarithm A^^^ has also good approximation for small values 
of \(\ for the same C = {e—z)/e: 



(39) 



A(2)(j 



3 



16 



with the same coefficients c, as in the case of A'-^K (The only difference is the 
opposite sign in the argument of the logarithm.) For large values, an extension to 
a wide range in the complex plane can be similarly realized with 



log(z) 



1 



(40) 

The resulting function 

(41) Asiz) = A(2)(z) - A(2)(3) sa A'-^^z) + 20.0563555297533789 

is plotted in the right bottom part of figure [3j Function A3{z) is not periodic, and 
has cut from the branch-point z = e to the negative direction of the real axis. In 
the positive direction of the real axis, it grows to infinity, and this grow is very slow 
(slower than any finite combination of logarithms). 

The asymptotic representation for the Abel functions Ai and A^ can be inverted, 
using various combinations of A—e and ln(±(A— e)) as a small parameter. Different 
small parameter allows different representations for the super-exponentials Fi and 
F; it seems many of them give comparable speed and comparable prevision; at 
least they do not add much errors to the rounding errors at the complex < doube > 
implementation. The asymptotics of the previous section seems to be the fastest, 
although the careful comparison of efficiency of various asymptotic formulas may 
be subject for the future investigation. 

For the plotting of figure [3] the algorithms for Ai and A3 were implemented in 
C-|— I- with complex < doube > arithmetics. In order to verify the consistency of 
these algorithms to those for Fi and F3, the following agreements are considered: 



(42) 



D 



lAF 



(z) = Ig 



(43) i?3AF(^)-lg 



Ai(fi(z)) + z 
Ai(Fi(z))-z 

A3(F3(z)) + Z 



A3(F3(z))-Z 



Difa{z) = Ig 



D3Fa{z) = Ig 



Fi(^i(z)) + z 
Fi{Ai{z))-z 

F^iAsiz)) + z 



F3{A3{z))-z 



Two of them, namely, D = D\p^-p{z) and D = D-j,FK{,z) are shown in figure [4] with 
contours D^const. As for DyppJ\z) and I?3af(-z), they remain of order of 14 in the 
whole range of such a picture, so, they are not presented here. 
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Figure 4. The agreements D = Dipa{z) and D = D^Apiz) and 
in the complex plane z = x+iy; level Z) = l is shown with thick lines; 
the integer values are shown with thin lines; levels Z? = 2 is shown 
with black thick lines, levels D = 1,2, 10, 12, 14 are seen; symbols 
"15" and "11" indicate the ranges where Z5>14 and 10 > 13 > 12. 
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In figure |4j symbol "15" indicates the region where D > 14, and symbol "11" 
indicates the ranges where 10 < D < 12. 

The top picture in figure |4] shows, that in the central part, the implementations 
of functions Fi and Ai = F^^ are consistent within at least 14 decimal digits. In 
the right hand side, the branches of functions Fi and Ai do not match, and the 
agreement is poor. 

The good agreement indicates, that the algorithms above work close to the best 
precision achievable with the coinplex<double> variables. However, the range of 
validity of relation F{F~^{z)) — z is limited by the cut lines of the functions F and 
F~^. In such a way, the numerical tests confirm the efficient C++ implementation 
of the super-exponentials Fi and F^ and the corresponding Abel-exponentials Ai 
and ^3. 



5. Non-integer iteration 

Each of the pairs {Fi,Ai) and (^3,^3) can be used to construct the regular 
iteration of the exponential to base & = exp(l/e): 

(44) expg(z) = Fi(c + ^i(z)) 

(45) e^p^^,{z) = F3{c + A^{z)) 

These functions are shown in figure [5] for c= 1/2. For comparison, the function 
y = expij{z) is plotted with a thin curve. Visually, the thick solid curve looks like 
a continuation of the dashed curve. However, analytic continuation is not possible 
because z = e is a branch point oiAi{z). Similar visual effects are discussed in [5] for 
the case b = ^/2. In order to see that exp['^'j^(z) and exp['^g(z), at least for 6 = exp(l/e) 
and c=l/2, are pretty different functions, the difference of the "continuations" of 
these two functions, id est, 

(46) dqi3{x) = exp[^(^'(a:+io) - exp[,^^^' (x+io) 

in vicinity of a; = c is shown at the bottom picture of figure |6] 

The complex map of the two square roots of the exponential to base exp(l/e) are 
potted in figure[6]in the same notations, as in figures [2] and |3] The additional levels 
p = 0.692907175521155 and q = ±8.53 are plotted in order to reveal the behavior 
of the exp[°j^^'(z) in the strips along the cut lines ^{z) > 0.7, 5(z) ~ ±8.5 . The 
comparison of the pictures in figure |6j show how different the functions exp^*^]^^' and 
expjj'g^' are when considered away from the real axis. For the functions exp[^^^' (z) 
and expjj^g^' (z), in wide ranges of z, the relations 

(47) expffl(expffl(z))=z 

(48) expf3/^](expf3/^l(z))=z 

hold. As in the case of the square root of the logistic operator [TT], the ranges 
of validity of these equations do not cover the whole complex plane, and they are 
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Figure 5. Behavior of the two square roots (c = of the expo- 
nential to base 5= exp(l/e). Top picture: Dashed: ?/ = exp[^(^' (x^ 

Thin sohd: y - 



by (|44|); Thick sohd: y = exp[^^^'(a;) by (|45| 
expj(x) = . Bottom picture: y ~ ^{dqi3{x)), thick hne, and 
2/ = 9((iqi3(a;)), thin hne, in vicinity of a; = e by equation (46 1. 



different. In order to show these ranges, the agreements 

+ expb(z) 



(49) 



(50) 



^q3(2)=lg 



[1/2] / [1/2]. s 

exp^ [expl[ '(z) 



[1/2] / [1/2]. N 

exp^ [expl[ '(z) 



[1/2] 



(expiy'(z) 
(expS/^l(.) 



expb(z) 

-expb(z) 
- expf,(z) 



are shown in figure [rj In particular, exp\^^'^\x) can be considered as "true" square 

root of the exponential to base b = exp(l/e) for a; < e, while exp[^3^^(a:) can be 
considered at the "true" root at a;>e; but these square roots can not be combined 
into the same holomorphic function. 

The fractional iterations by (44 1, (44) can be evaluated for complex values of z 
and even for complex values of c; but only at integer values of c these two functions 
can be considered as holomorphic extensions of each other. The fractional iteration 
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-8 -6 -4 -2 2 4 6 8 10 12 14 IG 18 20 22 24 26 28 

Figure 6. Two different "square roots" of the exponential to base 



exp(l/e) by equations (44 1 and (45) in the same notations as in 
figures [2] and [3j 
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Figure 7. The agreements Dqi{z) and Dq^lz) by equations (491 
and (50) in the complex plane z — x+iy. 
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Figure 8. Tetrationals to base 6=10, & = e«2.71, 6 = 2, 6=1.5, 
6 = exp(l/e)ssl.44 (the same curve as Fi in figure[lj, and 6 = a/2~ 
1.41 versus real argument. 

provides a smooth (holomorphic) transition from the exponential at c = 1 to the 
logarithm at c = — 1, passing through the "square root" of the exponential at 2: = 1/2, 
the identity function at c = and the "square root" of the logarithm at c=— 1/2 . 
In a similar way, the complex iterations of a function can be considered. 

The non- integer iteration of the exponentials provides a set of functions that grow 
up faster than any polynomial but slower than any exponential. Such functions 
may find applications in various areas of physics and technology. Similar non- 
integer iteration for other functions (including the exponentials of different bases 
and factorial) were discussed recently [71 ^ (TUl [H] , but the peculiarity of the fixed 
points of the exponential at the base 6 = exp(l/e) required the special consideration 
above. 

6. Comparison of the tetrationals to different bases 

In this section, the tetrational Fi to base exp(l/e) is compared to tetrationals 
to various bases. In fig [8] the tetrational tetf, versus real argument is shown for 
6= 10, 6 = e« 2.71, 6 = 2, 6= 1.5, 6 = exp(l/e) w 1.44 and 6= \/2 sa 1.41 . The 
functions for 6>exp(l/e) are evaluated using the Cauchy algorithm described in 
[7]. For 6<exp(l/e), the regular iteration described in [3] is used. For 6 = exp(l/e), 
the tetrational is just Fi shown also in figure IT]). 
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Figure 9. Tetrationals tet;,(a;+i?;) = p+iq for various b, from left 
to right: 6 = 1.5, 6 = exp(l/e)~ 1.44 (the central part of map of Fi 
at the top picture in figures, and 6=\/2~l-41 in the x,y plane. 



At moderate values of argument or order of unity or smaller, the curves for 
b = 1.5, 6 = exp(l/e) and b = \/2 are very close. In order to see the difference, the 
complex maps are shown in figure|9]for b= 1.5, left, for 6 = exp(l/e), central, and for 
b—\/2, right. The central picture of figurejojis just a zoom-in from the central part 
of the top picture in figure [2] The efficient algorithms of the computation allow to 
plot all the figures with some reserve of resolution; at the online version, they still 
can be zoomed-in. 

Figure [9] indicates no qualitative change of the tetrational at small variation 
of the base in vicinity of value b = exp(l/e). In particular, within the loop p = 
3fJ(tetfc(z)) — 1 (this loop goes through the origin of coordinates in all the three 
pictures) , the zooming- in of the central parts of pictures in figure |9] is necessary to 
see the difference. 



7. Comparison of the non-integer iterates of the exponential to 

different bases 



10 



versus x for b 



The non-integer iterates expj^ [x) are shown in figure 
b = exp(l/e), b = \f2 and various values of c. At 1 < < exp(l/e) there are two 
c-iterates of the exponential. The one valid for x below the upper fixed point is 
shown with a dashed line; the one valid for x above the lower fixed point is drawn 
with a solid line. The curves for 6 = exp(l/e), c= 1/2 and c = 1 are the same as 
in top picture of figure [5j For the evaluation of the iterated exponential to base 
6 = e, the approximation by [8] is used (although the direct method by [7^ could be 
used instead) ; the two iterated exponentials to base b = \pl were evaluated with the 
algorithms described in i9^. In the last case, in the range 2<a;<4 both the iterated 
exponentials are valid, and the deviation between the two iterated exponentials is 
of order of 10~^^; however, even in this case they are not holomorphic extensions 
of each other. (Being plotted on a paper, the distance between the two curves for 
the same c is small not only in comparison to the size of an atom, but also small 
being compared to the size of atomic nuclei.) 

8. Conclusion and prospects 



We present some theory of regular iteration and apply it to the case /(z) = e^/''. 
We extract a quite efficient algorithm to calculate the two super-logarithms to base 
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6"'^/'' compared with various other (standard and separate) methods. We suggest an 
efficient new non-polynomial approximation for the two super-exponentials to base 
ei/<=. 

One of the two super-exponentials, the tetrational Fi = tetoxp(i/o)i is holomor- 
phic in the range C \ e M : a; < —2} and strictly increasing along the real axis 
> —2. It asymptotically approaches the limiting value e. The function approaches 
the same value e also in any other directions, i.e. at large values of \z\. The jump 
at the cut z<—2 reduces to zero, as \z\ — > co. 

The other super-exponential -F3 is an entire function. Along the real axis it is 
strictly increasing from the limiting value e at —00 to infinity, growing faster than 
any exponential. Outside the positive part of the real axis, ^3(2:) approaches e at 
\z\ —^00 in a similar way as Fi does. 

Efficient calculation algorithms and portraits for bases b > e^/° and 1 < b < 
e^/° were already provided in [3 [3]; so the whole range 6> 1 is now covered. The 
plots of the tetrational tet;, for 6 = e-'^/°« 1.44 look similar to those for b— 1.5 and 
those for b — 1.41 ; one may expect that at any fixed value of z from some 

range, the tetrational tetf,(z) is a continuous function of b at least for 6> 1. This 
raises the following question: 

Are the holomorphic tetrationals constructed in 9 , here and in ^17j (generaliza- 
tion of 6J which is conjectured to be the super-exponential in [7]) — which together 
cover the base range (l,c») — analytic as a function of the base 6, particularly in 
the point b = e^/°? If so, what is the range of holomorphism? If not: can one 
obtain an operation (6, z) 1— > tetf,(z) defined for 5 in a vicinity of e^/° such that for 
each b the function z 1— >■ teti,{z) is a real-analytic tetrational on (—2, 00) and the 
function b 1— tetfc(z) is holomorphic for each z? 

There is a similar bifurcation base b « 1.6353 for the tetrational tet{,{z) as the 
bifurcation base e^^'^ is for the exponential (i.e. where the two fixed points change 
into no fixed point). One could apply the same methods we used to obtain the 
super-exponentials to also obtain a super-tetrationals/pentationals. 
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